Superfluid 3 He in a restricted geometry with a perpendicular magnetic field 
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We theoretically investigate the role of surface Andreev bound states (SABSs) on the phase 
diagram and spin susceptibilities of superfluid 3 He confined to a restricted geometry. We first 
explicitly derive the dispersion of the SABS in 3 He-B in the presence of a magnetic field, where the 
Majorana Ising spin and the spin susceptibility contributed from the SABS are associated with the 
5*0(3) order parameter manifold. Subsequently, based on the quasiclassical Eilenberger theory with 
Fermi liquid corrections, we discuss the nonlinear effect of a magnetic field on the SABS, where the 
magnetic field is perpendicular to the specular surface. It is directly demonstrated that a gapped 
SABS strongly enhances the magnetization density and spin susceptibility at the surface, compared 
with that in the normal 3 He. To capture the characteristics of the SABS, we show the field- and 
temperature-dependences of the spatially averaged susceptibility which is detectable through NMR 
experiments. It turns out that the contribution of the SABS leads to nonmonotonic temperature- 
dependence of the spin susceptibility. Furthermore, we present the superfluid phase diagram, where 
the B-phase undergoes a first-order (second-order) phase transition to A-phase or planar phase at 
low (high) temperatures. 

PACS numbers: 67.30.H-, 67.30.ht, 67.30.ef, 74.20.Rp 
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I. INTRODUCTION 

The study of spin-triplet p-wavc superfluid 3 Hc in a re- 
stricted geometry has a long history since the mid 1970s. 
It is known that the most symmetric superfluid phase, 
called the B-phase, is a ground state in the bulk 3 He 
at low temperatures, because the isotropic energy gap 
gains the more condensation energy than other compet- 
itive phases*^— In a restricted geometry, however, sur- 
faces give rise to the pair breaking, which squashes the 
isotropic B-phase clliptically. As two surfaces get close to 
each other, the B-phase at last undergoes a change to the 
A-phase which has point nodes at the north and south 
poles of the Fermi sphere. The A-B phase transition in- 
duced by a surface boundary condition was directly ob- 
served in NMR experiments and a torsional oscillator^— 
where superfluid 3 Hc is confined to a geometry with sub- 
micron thickness. 

Theoretical studies for understanding the pair breaking 
effect and the n-texture on surfaces were initiated by the 
analysis based on the Ginzburg-Landau theory^ - — Be- 
yond the theory which does not take account of the infor- 
mation on quasiparticles, the quasiclassical Eilenberger 
theory provides a tractable and quantitative scheme to 
study the interplay of the pair potential and quasiparticle 
statesi 25 ' 26 It was found by Buchholtz and Zwicknagl^I 
in 1981 that from the microscopic point of view, midgap 
bound states emerge on a specular surface of the B-phase. 
The midgap state was more explicitly discussed in Ref.l28l 
which analyzed a p-wave polar state as an exactly soluble 
model. Based on the quasiclassical theory, the quantita- 
tive phase diagram and the finite size effect of the su- 
perfluid 3 He in a restricted geometry have been clarified 
in Refs. I29U32I which underline the role of the SABS on 
thermodynamics. Note that these previous works have 
not taken account of the effect of a Zeeman magnetic 



field. Several experiments have observed the pair break- 
ing effect and the enhancement of the surface density of 
states due to the midgap states 3 - - — concurrently with 
theoretical works. It is worth mentioning that at the 
present time, the midgap state is recognized as a family 
of the Andreev bound state, called the surface Andreev 
bound state (SABS)r^. Since this midgap bound state 
exists only if the superconducting pair potential changes 
its sign^ it ubiquitously appears in various physics 
systems, such as superconducting junctions^ super- 
conducting vorticeS ) ' unconventional superconduc- 
tors^^— and Fulde-Ferrell-Larkin-Ovchinnikov super- 
conductors In addition, the same physics is shared 
with the solitons which emerge in polyacetylene^i - — the 
incommensurate spin-density wave ; 55 ' 56 and the stripe 
state in high T c cuprates^i 

Recently, the study of the SABS in superfluid 3 He- 
B has been rekindled by the understanding of the di- 
rect correspondence to a bulk topological invariant. The 
chiral symmetry constructed from the time-reversal and 
particle-hole operations allows us to introduce the three- 
dimensional winding number which behaves as a topolog- 
ical invariant in 3 He-B^2r— Hence, the gapless SABS in 
3 Hc-B is protected by the topological invariant associated 
with the time-reversal symmetry. Furthermore, the chiral 
symmetry ensures that the SABS behaves as a Majorana 
fermion which is a particle equivalent to its own anti- 
particlei^— As a consequence of a Majorana fermion 
and the topological protected bound state, it is unveiled 
that the multifaceted SABS is insensitive to the density 
fluctuation and exhibits Ising-like isotropy of magnetic 
response, that is, what is known as the Majorana Ising 
spin i 69 ! 71 ' 72 Recently, the gapless cone spectrum has been 
observed by measuring the transverse acoustic impedance 
on surface of 3 He-B with controlled specularity . 73 ' 74 

The Majorana fermion and Ising spin may survive even 
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when a magnetic field is applied along a direction parallel 
to the specular surface, where the time-reversal symme- 
try is broken. In a weak magnetic field comparable with 
the dipolar field, the B-phase stays topological as a con- 
sequence of the protection with a hidden Z 2 symmetry, 
that is, called the symmetry protected topological or- 
der^ This ensures the chiral symmetry and protects the 
gapless Majorana cone even in the presence of a magnetic 
field. In the case of a perpendicular field, however, an in- 
finitesimal magnetic field is always able to open a finite 
energy gap in the surface Majorana cone. Within the 
linear regime of the applied field, the spin susceptibility 
with respect to a perpendicular magnetic field is consider- 
ably enhanced, compared with that in the normal 3 He.— 
This enhancement is expected to make the SABS de- 
tectable in NMR experiments. Nevertheless, there have 
never been any studies which quantitatively and micro- 
scopically discuss thermodynamics and spin susceptibil- 
ity in the nonlinear regime of the magnetic field. Note 
that in the magnetic field regime much stronger than the 
dipolar field (^30G), the Majorana Ising anisotropy dis- 
appears and the thermodynamics is independent of the 
orientation of the applied magnetic field j 70 ' 75 ' 76 

Hence, the purpose of this paper is to clarify the role 
of the SABS on thermodynamics and spin susceptibility 
in superfiuid 3 He restricted to a slab geometry. As dis- 
played in Fig. [TJ we set 3 He to be sandwiched by two 
specular walls and a magnetic field is applied along the 
z-axis which is normal to the wall. First, we clarify that 
the Majorana nature of the SABS is associated with the 
50(3) order parameter manifold of 3 Hc-B. Based on the 
quasiclassical Eilenberger theory taking account of Fermi 
liquid corrections, we quantitatively discuss the thermo- 
dynamics and the enhancement of the magnetization den- 
sity due to a gapped SABS in a slab geometry. We here 
present the quantitative superfiuid phase diagram in 3 Hc 
restricted to a slab, where it is emphasized that the Fermi 
liquid correction plays a critical role on determining the 
A-B phase boundary induced by the magnetic field. To 
capture the characteristics in experiments, we clarify the 
field- and temperature-dependences of the spatially av- 
eraged spin susceptibility, which is detectable in NMR 
experiments. In particular, we emphasize the nonlinear 
effect of the applied magnetic field. It turns out that the 
spatially averaged spin susceptibility in a 3 He with a sub- 
micron thickness exhibits the non-monotonic behavior on 
the temperature. 

In the following section, we explicitly derive the disper- 
sion of the surface bound states in the B-phase parame- 
terized with the n-vector and ip which are the parameter 
of the 5*0(3) order parameter manifold. Here, we dis- 
cuss the relation between the £0(3) manifold and spin 
susceptibility contributed from the surface bound state. 
Then, we move on to the quantitative calculation based 
on the quasiclassical Eilenberger theory. In Sec. Ill, we 
describe the details of the quasiclassical formulation for 
superfiuid 3 Hc in a restricted geometry, where the self- 
consistent framework takes account of the Fermi liquid 



corrections. The self-consistent solutions for the pair po- 
tential, magnetization density, and local density of states 
are presented in Sec. IV, where the contribution of the 
surface bound state to magnetization density is under- 
lined. The complete phase diagram is proposed in Sec. V, 
where we demonstrate that the phase boundary is sensi- 
tive to the Fermi liquid corrections. We also present the 
non-linear effect of the Zeeman magnetic field on the spin 
susceptibility. The final section is devoted to conclusion 
and discussion. The details on the derivation of the dis- 
persion of the SABS are described in Appendix A. The 
numerical procedure and boundary condition are given in 
Appendix B and Appendix C shows that the quasiclassi- 
cal Eilenberger equation for 3 He-B is invariant under an 
SO (2) rotation in spin and orbital spaces, which short- 
ens computation time. Throughout this paper, we set 
h = feu = 1 and repeated Greek (Roman) indices imply 
the sum over x, y, z (y, \). 



II. SURFACE BOUND STATE AND 
MAJORANA ISING SPIN IN THE B-PHASE 

A. The B-phase order parameter 

Let us start with the mean-field Hamiltonian density 
for superfiuid 3 He with the mass M in the 4x4 Nambu 
representation, 

Kin, r 2 ) = 5{r 12 ) [e(ri)r, + V] + E(n, r 2 ). (1) 

The single-particle Hamiltonian density is e(r) = 
— V 2 /2A7 — Ep and the Zeeman energy V_ is given by 

_y, (2) 

with the Fermi energy E-p = k F /2M, and magnetic mo- 
ment of 3 He nuclei /J, n . We also introduce the Pauli ma- 
trices Up (t^) in the spin (particle- hole) spaces. The 
self-energy matrix S(ri,T"2) consists of the Fermi liq- 
uid correction E FL and the pair potential A, that is, 



z 




D 



FIG. 1: (Color online) Schematic picture of the system which 
is considered in this work. The superfiuid 3 He is sandwiched 
by two specular walls whose distance is set to be D. In and 
after Sec. Ill, a magnetic field is applied to the z-axis which 
is normal to the wall. 
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S = E FL + A, and 



A(n,r 2 ) 



A(n,r 3 ) 
-A*(n,r 2 ) 



(3) 



where A is a 2 x 2 matrix in the spin space and the 
spin triplet and p-wave symmetries require A a b(ri, r 2 ) = 
A ba (r 1 ,r 2 ) and A ab (r 1 ,r 2 ) = -A ab (r 2 ,r 1 ). 

Without loss of generality, a pair potential for a 
spin-triplet supcrfluid is expressed with the d-vector as 
A(fe,r) = / dr 12 e-* kr ^A( ri ,r 2 ) = ia^yd^r). The 
supcrnuid 3 Hc-B phase is known to be most symmet- 
ric among possible order parameters, which is invariant 
under the joint rotations of three-dimensional spin and 
orbital spaces, SO(3)l+s- The order parameter is the 
eigenstate of the angular momentum operator composed 
of the spin and orbital angular momentum, S and L, 



= + S v R Uil {n, if), 



(4) 



implying the spontaneously broken spin-orbit symme- 
try, SO(3)l-s, in addition to the ordinary U(l) gauge 
Here, i? Mi ,(n, (p) describes the relative rotation ma- 
trix between spin and orbital spaces originated from the 
SO(3)l~s manifold, where n and ip denote the rota- 
tion axis and angle. Then, the general form of the or- 
der parameter of the supcrnuid 3 He-B is described as 



d li {k,r) = d tlv {r)k v 



d^ir) 



jo 



R,j, v (n,<p)A„(r) 



(5) 



This order parameter also describes the planar phase 
with A z = and squashed B (or B 2 ) phase with A z ^ 
A x = A y Jul Since the planar phase is energetically de- 
generate with the A- phase at the weak coupling limit, 
the order parameter in Eq. © takes account of all possi- 
ble phases stabilized in the presence of a magnetic field. 

Here, to derive the SABS from the general form in 
Eq. ([5]), let us consider the situation where a single spec- 
ular surface is set at z = 0. The energy eigenstates of 
the Hamiltonian in Eq. ([TJ are obtained by solving the 
following eigenvalue equation, 



J dr 2 n(ri,r 2 )ip(r 2 ) = E(p(r 2 ), 



(0) 



which is called the Bogoliubov-de Gennes (BdG) equa- 
tion. Here, the wavefunctions obey the normalization 
condition, J drip^ (r)cp(r) = 1. To solve Eq. analyti- 
cally, in this section, we ignore the Fermi liquid correc- 
tion, that is, I] pl = 0. The effect is discussed in the subse- 
quent sections with the quasiclassical Eilcnberger theory. 
In addition, the pair potential A is assumed to be spa- 
tially uniform A(fc, r) = A(k). It is often convenient to 
utilize the alternative description with U(n,tp) <S SU(2) 
as U{h 1 ip)a v U\h 1 ip) = a^R^(h,ip). Using this SU(2) 
matrix, the B-phase order parameter in Eq. ([5]) reduces 
to 



where U T denotes the transpose of a matrix U and 
Ao(fc,r) = ie^o^Oy A AI (r)fc M is the simplest expression 
of the B-phase order parameter. 

Using the Andreev approximation which holds within 
the weak coupling regime kp£ = 2Ep/A 3> 1, the BdG 
equation ([6]) reduces to the Andreev equation 

[-iav F cos 9 k d z T z + V + A(fc F , Q )] <£a(z) = E<p a (z), (8) 

where <f±(z) describes the slowly varying part of 
quasiparticle wavefunction <p(r), that is, ¥>(r) = 
X^o^h C a ip a (z)e lkF ' a ' r with the Fermi velocity vp and 
kp.a = fcp(cos0fc sin6>fc, sin^fc sin6>fc, a cos k 9 k ). 



B. Majorana Ising spin 

First, we solve the Andreev equation ([8]) in the absence 
of a magnetic field, H — 0, where the resulting equation 
becomes equivalent to that for spinlcss chiral p-wave su- 
perconductors^ As described in Appendix A, the bound 
state solution with |i?(fcii)| < Ao has the energy disper- 
sion linear on the momentum fen = {k Xl k y ) as 



£ (fc||) = ±^f|fc|| 



(9) 



This expression is independent of the orientation of n 
and the angle tp. The corresponding wavefunctions for 
the quasiparticles bound at at z = are given by 

vS„( r) = N k e ik rnf(k uZ ) U {n, V )^ ± {cj> k ), (10) 

where N k is the normalization constant and hi = 
diag(C7, U*). In Eq. ([TO]), we also set f{k±,z) = 
sin (k±z) e~ z /^ with k± = ^/fc| — k^ and 



*±(0fc) 























e 2 





















(11) 



The wavefunction <p^^ corresponds to the positive en- 
ergy solution of -Eb(fc|| ) and ip^ ^ is the negative branch. 

The particle-hole symmetry T x H(k,r)T_l, = — T-L*(— k,r) 
ensures the one-to-one correspondence between the two 



branches as ip^ ^ 



(+)* 



(* r ,\E' i ,«'|,4'{) T in spin- 



A(fc,r) = U(n,cp)A (k,r)U T {n,<p), 



(7) 



The quantized field \I/ 
triplet supcrfluids can be expanded in terms of the pos- 
itive energy states of the SABS with E(k^) > and 
cp kll (r) in addition to continuum states. For low tem- 
perature regimes T<cA , the field operator can be con- 
structed from the contributions of only the SABS as 

*(»•) ~ Efc„ [VoSi ( r )%ii + r x f^ k *(r)vi^ where 7? fe|| 
and rj\, denote the Bogoliubov quasiparticle operators. 
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Then, the field operator contributed from the SABS 
obeys the self-conjugate Majorana condition, 



*M (r)RJ _[*M (r) ]t 



(12) 



where ^^(r) = U'J b (n,(p)^' b (r). Using the wave- 
function of the SABS, *«(»-) is given as *'(r) = 



D ife|| -T-y +i<p h /2 



h.c.]Z/(n,</j)*' fe . Here, the 



spin quantization axis is changed from the z to x- 
axis, where W, change to St', W, 3?^. = 

[cos^-, sin^-, — cos^, — sin^ L ] T . 

Once Eq. (fT2"j) holds, it is straightforward to prove 
that the Majorana fields behave as the Clifford 
algebra, {^(n), tfg*(r a )} = 25 0>6 ^(r 12 ). For the 
case of n || £, the local spin operator 5 M (r) = 
| x l , 3 I (r)(o' M ) a ; ) \E'f ) (r) with the Clifford algebra results in 
the Ising-like anisotropic form as S = (0,0, S^ 1 ) = S M , 
where S^ 1 = -~^l^(a^) ab ^. For an arbitrary configu- 
ration of (n, </?), the local spin operator contributed from 
the SABS results in 



S fl (r)=R )iz (n,cp)S™(r) 



(13) 



spectrum in Eq. (|9|), where D denotes the thickness of 
the sample. Then, the hybridization of the two SABSs 
exponentially splits the zero energy state at |feii| = 
with quantum oscillation on the scale of fcp 1 as 6E(k\\ = 
0) ~ e~ D ^ sin(kFD)^— In the quasiclassical Eilcn- 
berger theory, the quantum oscillation term vanishes^ 
In addition to the splitting due to the quasiparticle tun- 
neling, the finite size of the system with the thickness 
D = 0{^) gives rise to the pair breaking effect, which may 
also stimulate the deviation of the gaplcss spectrum. The 
distortion of the gapless Majorana cone due to the quasi- 
particle tunneling and pair breaking effect may break 
the Majorana Ising nature of the surface bound states. 
The numerical analysis on this issue will be discussed in 
Sec. IV. 

Now let us turn to the case of a finite magnetic field 
Hy^O. As described in Appendix A, the dispersion of the 
SABS is given as 



E(ku)=±\\E (k 



H a H£ z (n 1 (p) 



(16) 



The direction of the Majorana Ising spin reflects 
the 50(3) order parameter manifold (n,<p). Using 
Eq. (| 13[> . the dynamical spin susceptibility becomes 
X^(ri,r 2 ;uj) = xf z {ri,r 2 -,uj)R^ z (h,(p)R vz (n,ip). This 
implies that magnetization and susceptibility originate 

from Majorana Ising spins Sjf(r) and Xzz\ r ii r 2',w) = 
(Sf (ri)Sf (r 2 )) u through the 50(3) matrix R^(n,ip). 
The property of Xzz { r ii r 2] was discussed in Ref. HI, 

M 

To understand the orientation of the Majorana Ising 
spin in Eq. (|13[) , it is convenient to introduce the ^-vector 
in Eq. ([TBI , the dcfinitio n 70 i 75 of which is 



l M (n, <p) = h v R vli {n, ip) 



(14) 



The orientation of an applied magnetic field is denoted 
by h v — Hy/H. Then, it turns out that the £ z (n,ip) 
describes the projection of the Majorana Ising spin S(r) 
in Eq. (|13[) onto the orientation of the applied magnetic 
field H as 



and the wave functions are obtained from Eq. (|A13[) with 
a±(fc||) = 



|(1 ± 



|Bo(fc||)| 



j I). The resulting dispersion in 
Eq. (p~6|) implies that the energy gap of the surface state 
depends on the ^-vector as 



(17) 



min £?(fc||) = [i n H £ z (n,ip) . 



For h = z, since l z — h z , the dispersion in Eq. (flBT) is 
consistent with the previous works in Refs. I?32ll69ll7lll72l . 
Equation (fT6|) indicates that for £ z ^ 0, the Majorana 
Ising nature of the SABS disappears and an arbitrary ori- 
entation of the magnetic field opens a finite energy gapi 7 -^ 
The Ising anisotropy is also consistent with Eq. (fl"5|) de- 
scribing the relation between the orientation of the Ma- 
jorana Ising spin and the applied field. 



h ■ S(r) 



(15) 



Figure [2] depicts the schematic picture for S, H, and £ z . 
For £ z = 0, the Majorana Ising spin S is perpendicular to 
the applied magnetic field, which implies that the SABS 
does not contribute to the magnetic response. However, 
the SABS may be responsible to H when £ z ^0. 

The gapless spectrum of the SABS is protected by the 
nontrivial topological invariant defined in the bulk re- 
gion of the B-phase i 60 ' 61 i 70 As two specular surfaces at 
z = and D get close to each other, however, the in- 
terference between the SABSs distorts the surface cone 




FIG. 2: (Color online) Schematic picture on the relation be- 
tween i z and S, where S and S M denote the orientation of 
the Majorana Ising spins for an arbitrary n and S for n || z. 
The z-axis is normal to the specular surfaces as displayed in 

Fig.m 
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Effect of a Zeeman magnetic field and the 
dipole interaction 



The quasiclassical Green's function g for spin-triplet su- 
pcrfiuids is described in the particle-hole space as 



It is important to mention the relation between the 
energy gap of the SABS in Eq. (jTT)) and the stable con- 
figuration of (n,ip). The magnetic field energy density 
within the Ginzburg-Landau theory is given as 

2 2 
Fr = -es [H^R^{n,Lp)s v ] = -e s H 2 i z (h,ip) , (18) 

where es = £,o(xn — Xb) > 0.— The minimization condi- 
tion of -Fh, that is, \t z (n, tp)\ =1, opens the maximum en- 
ergy gap of the surface Andreev bound state in Eq. (fTT|) , 
min \E(k\\)\=(M u H. 

On the other hand, the dipole energy favors to align n 
to the i-axis. Hence, it turns out that a Zeeman mag- 
netic field perpendicular to the surface always open a 
finite energy gap in the SABS, min |i?(fe||)| = fi n H, be- 
cause both the magnetic field energy and dipole energy 
favor h || z. In the case of a parallel magnetic field, 
however, the dipole interaction energy may be competi- 
tive to the magnetic field energy. Actually, it is demon- 
strated in Ref. [lO that the gapless SABS with £ z may 
be protected by a hidden Z 2 symmetry which is pre- 
served under a magnetic field weaker than the dipolar 
field is regarded as the symmetry protected topological 
phase. Here, l z behaves as the symmetry protected topo- 
logical order. However, the topological phase transition 
with the spontaneous symmetry breaking takes place at a 
magnetic field comparable with the dipolar field beyond 
which l z = \ is realized to minimize the Zeeman magnetic 
energy. Since the magnetic response becomes isotropic in 
the high field regime, without loss of generality, the fol- 
lowing section will focus on the simple situation where 
the magnetic field is applied along the surface normal. 

Note that the angle (p is locked by minimizing the 
dipole interaction to 17 ' 83 " — 



l (A ± (z)A||(z)) 
4 (Af(z)) 



(19) 



A, 



A y and A 



A, and 



where we set Ay 

T> Jif ' ' ' d z - Equation (|T9| depends on the ratio of the 
pair potentials which are distorted by the nonlinear effect 
of the Zeeman magnetic field. 



III. QUASICLASSICAL THEORY FOR 
SUPERFLUID 3 HE 



(21) 



where f = f (k,r;iu} n ) and fp = fp(k,r;iu} n ) denote the 
spin-singlet and -triplet components of the quasiclassical 
Green's function, respectively. Note that g satisfies the 
normalization condition g 2 = —tt 2 t q . Here, we introduce 
<7o and Tq as the unit matrix in spin and particle-hole 
spaces. 

The evolution of the Nambu-Gor'kov Green's func- 
tions G is governed by the Nambu-Gor'kov equation, 
[— icu n + 2£\G = tq. Following the standard procedure, 
the quasiclassical Green's functions g(k, r; ie m ) obeys the 
so-called Eilcnbcrger equation,— 



iuj„T z - <S(fe, r) - v, g(k, r; iw n ) 

+ivp(k) ■ V<?(fe, r; iui n ) = 0. 



(22) 



The quasiclassical Green's functions must satisfy a con- 
straint given by the normalization condition, g 2 = —tt 2 t . 

The Fermi velocity vp is given as v-p(k) = vpk on the 
three-dimensional Fermi sphere. The Zeeman magnetic 
field is included in Eq. (f2"2")l as 



1 



1 



rp a, — z 



1 



1 + F} 



a 







(23) 



where Fq is one of the Landau parameters which de- 
scribes the enhancement of the spin susceptibility, as 
mentioned below. The 4x4 matrix S_ describes the qua- 
siclassical self-energies obtained from S_(k,r) w AE(fc = 
k F k,r)r z , where E(fc,r) = / dr 12 e tk - ri2 'S{r 1 , r 2 ). The 
quasiclassical self-energy matrix consists of the Fermi liq- 
uid correction in the diagonal elements and the pair po- 
tential d,,, 



S(k,r) 



ia y a,j,d* 



V >1 



(24) 



where we set vq = Vo(k, r), v„ = Vu(k, r), and = 
dfj,(k,r). We also introduce the notation, v'^ „ = 

The Fermi liquid corrections vq and are associated 
with the quasiclassical Green's functions go and g v as 



(fc, r) = Y, A i ( p e(k ■ k')g (k', r; iu) n ) 



,(25a) 



The quasiclassical Green's function g is obtained from 
the Nambu-Gor'kov Green's functions G with the Mat- 
subara frequency tu n = (2n + 1)ttT (n€Z) and the quasi- 
particle renormalization factor A as 



g(k,r;iuj n ) 



dt; k T z G(k,r;iLo n ). (20) 



-E a 



r) = V4 a) (P e (k ■ k')g^k',r;iLu n ))„ ,(25b) 



where {'")u denotes the Fermi surface average and 
Matsubara sum: (• • • ) k n =T£ K | <Eo / g. The Fermi 
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liquid corrections are expanded in terms of the Lcgcn 

(s) fa) 

dre polynomials Pi. The coefficients A\ and A e an 
the symmetric and antisymmetric quasiparticle scatter 
ing amplitudes, which are parametrized with the Lan- 
dau's Fermi liquid parameters,— -F?' a , through 



AT = 



1 v 



i + F; ,a /(2£ + i) 



(26) 



The I = (£ = 1) channel of the symmetric part in the 
Fermi liquid corrections couples to the density distribu- 
tion (mass current density J™(r)) and the antisymmetric 
part in the £ = and 1 channels arises from the magneti- 
zation density M M (r) and spin current density J^(r), 
respectively. They are defined with the quasiclassical 
Green's functions go and as 



M M (r) = M K 



H 



1 



pL n H 



J™(r) = 2v F N F (k„g (k,r 



\ I k' 



,(27a) 



(27b) 



(27c) 



where the magnetization in the normal state of 3 He is 
given by M N = 2/j,lN F H/(l + F^) = xnH and N F denotes 
the density of states at the Fermi energy in a normal 
Fermi gas. The spin current density J® describes the 
flow of the spin component along the f^-direction. As 
we will emphasize below, the parameter Fq coupled with 
the magnetization density strongly affects the qualitative 
feature of the A-B phase transition induced by a magnetic 
field. 

The pair potentials d tt (k,r) in spin-triplet supcrfluids 
are obtained from the gap equation with an attractive 
interparticle interaction V(k,k), 



df,(k,r) = (Vik^k^f^k.r^iujn) 



(28) 



where the pair interaction V is assumed to be invariant 
under the 5'0(3)i,x SO(3)s rotational symmetry in spin 
and orbital spaces. Hence, using the form V(k, k ) = 
igk^ fc' with the coupling constant g > and the form of 
the B-phase order parameter in Eq. ([5]), the gap equa- 
tions for ci M „(r) are 



(29) 



The coupling constant g is related to the transition tem- 
perature T c q in the bulk, which is given by the linearized 
gap equation at T~T c q as 



- = 7rT c0 V "j - 7 , 
9 , t-T 1 Wn.c 



(30) 



where ui c is the cutoff frequency and w„ iC denotes the 
Matsubara frequency at T~T c0 . 

In realistic situation of 3 Hc, the magnetic dipole in- 
teraction arises from the magnetic moment of 3 He nu- 
clei which reduces the SO(3)l x SO(3)s symmetry to 
SO(3)l+s- In the presence of a perpendicular magnetic 
field which we consider here, however, as discussed in 
Sec. II, the dipole interaction merely locks the angle ip 
to Eq. (fl"9)) and the contribution to the thermodynamics 
is negligible. Hence, for a perpendicular field, the order 
parameter reduces to 



(31) 



which corresponds to the case of n || z in Eq. ([5]). The 
effect of the dipole interaction becomes crucial in the case 
of a weak magnetic field parallel to the surface, which will 
be discussed elsewherei 70 i 86 

In summary, the Eilenberger equation (|2"2"j) coupled 
with Eqs. (|2"5|) and (|2T))) through the quasiclassical self- 
energies in Eq. (|24|) provides the closed form of the self- 
consistent equations for the quasiclassical Green's func- 
tions g and the meanfield potentials i/ Q , v^, and d^. In 
Appendix B, we describe in details the calculated sys- 
tems, boundary conditions, and the procedure for nu- 
merical calculations. 

As shown in Appendix C, the Eilenberger equation 
(f2"2l for 3 Hc-B is invariant under the SO(2)L^ + s r rota- 
tion, when the Zeeman magnetic field is applied along 
the surface normal. The symmetry leads to the one-to- 
one correspondence of the quasiclassical Green's function 
between two points k and O^k on the Fermi sphere, 
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(O^k, z; iu) n ) = U_-2g{k, z; iu n )U_ 2 



(32) 



where is an SO(2) rotation matrix about the z-axis 
and U 2 is the 4x4 matrix which describes an SU (2) rota- 
tion associated with O^- 2 ' (for the details, see Appendix 
C). This relation through the SO{2) rotation is useful for 
shorting the computation time of the selfconsistent cal- 
culation. Once we calculate g(k, z; iLu„) along the path 




FIG. 3: (Color online) One-to-one correspondence between 
two points on the three-dimensional Fermi sphere. The qua- 
siclassical Green's function at an arbitrary Ok belonging 
to the path (ii) is obtained by the SO(2)l z +s z rotation of 
g(k, z;iuj n ) calculated on the path (i). 
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(i) displayed in Fig. |3l the Green's function g for all fc is 

given by the symmetric relation in (|32|) with g(k, z; iuj n ). 

Throughout this paper, we use the set of the Fermi 
liquid parameters, F S = 9.3, F a = -0.695, f\ s = 5.39, 
and Ff = — 0.5ii The cutoff frequency on the Matsubara 
sum is taken to be uj c = 20ttT c q for low temperatures and 
1607rT c o for high temperatures. All length and energy 
scales arc in a unit of the coherence length in quasiclas- 
sical formalism, £o=Vf/ttT c o, and 7rT c o. 



IV. ORDER PARAMETERS, LOCAL SPIN 
SUSCEPTIBILITIES, AND SURFACE BOUND 
STATES 

A. Distortion of the B-phase order parameter and 
magnetization 

First of all, in Fig.^fla), we summarize the spatial pro- 
files of A||(z) and Aj_(z) for various thickness D where 
we fix H = and T = 0.2T c o. In the vicinity of the 
specular surface, the A z component is suppressed by the 
pair breaking effect and the parallel components remain 



isotropic, that is, A 3 



A || and A K 



A i . For 



a large D, e.g., D = 40£o, the isotropic B-phase order 
parameter with A|| = Aj_ appears around the middle re- 
gion z/D ~ 0.5. It is continuously turned to the pla- 
nar state with A z =0 at z = across the squashed B- 
phase with A|| > A z . As D decreases, the pair break- 
ing effect at the surface occurs even in the central re- 
gion, which elliptically squeezes the order parameters. 
The squashed B-phase undergoes a second-order phase 
transition to the planar or A-phase at the thickness 
D ps 9.6£ = D CI1 (H = 0)^2. when the magnetic field 
is absent. 

Figure 0Jb) shows the spatial profiles of Am (z) and 
A_l(» for /i n if/7rT c o = 0, 0.061, and 0.096 and L> = 20£ - 
It is seen that the magnetic field H || z as well as thick- 
ness D squeezes the B-phase order parameter elliptically, 
leading to A|| > Aj_. As we will mention in the subse- 
quent section, the squashed B-phase undergoes the first- 
order phase transition to the planar or A-phase at low 
temperature. 

The thickness-dependence of the local magnetization 
density M m (z)/Mn defined in Eq. (|27a[) is summarized in 
Fig. 03a). In the case of a large D/^q, the magnetization 
around the central region, e.g., z/D = 0.5, is strongly sup- 
pressed, compared with that in the normal 3 He. In the 
thermodynamic limit, D — >oo, the ratio of the magneti- 
zation between the B-phase and normal phase is obtained 
as X zz /xn = 2(1 + F a )/(3 + 2F a ) at T = 0^, which im- 
plies that Xzz~0-38xn for Fq =— 0.695. It is important 
to mention that for a large D/^q, a magnetic field per- 
pendicular to the surface enhances the low-temperature 
spin susceptibility on the surface, where M z (z = 0) > A/n 
and M z (z~D/2) <Mn- The enhancement is closely as- 
sociated with the energy spectrum of the surface bound 
states, which will be clarified in the subsequent subsec- 




FIG. 4: (a) Spatial profiles of A||(z) (solid lines) and Ai(z) 
(dashed lines) for D/^ = 10, 20, and 40 at H = and 
T = 0.2T c o. The horizontal axis is scaled with D. (b) Field 
dependence of Ay(z) (solid lines) and A±(z) (dashed lines) 
for D = 20£o and T = 0.2r c0 , where n n H /-kT cQ = Q, 0.061, and 
0.096. 



tion. As D approaches the critical value D cr i(0), the B- 
phase continuously changes to the planar phase through 
the squashed B-phase, where the spin susceptibility in 
the planar phase is indistinguishable from that in the 
normal state, Xz lanal =Xn- Hence, in this regime, the en- 
hancement ceases to exist and the magnetization density 
flattens due to the strong distortion of the B-phase order 
parameter. 

The local magnetization density feedbacks the effective 
magnetic field through the Fermi liquid corrections. This 
gives rise to a nonlinear effect of the Zeeman magnetic 
field. Since the distorted B-phase is not accompanied 
by the mass flow, the quasiclassical selfencrgics v(k, r) 
are composed of the local magnetization density M^(r) 
and the supcrfluid spin flow J s ^ v {r), which changes the 
Zeeman energy term to 



1+^ 



i^ tt (fc,rK,(33) 



where H° E (k,r) denotes the magnetic field deviated by 
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FIG. 5: (a) Spatial profiles of the local magnetization den- 
sity M m (z)/Mn, corresponding to the ratio of the local spin 
susceptibility Xzz( 2 )/xn, for various D's at T = 0.2T c o and 
fj, n H/TvT c o = 0.0122. (b) Field-dependence of Xzz( 2 ) on the 
surface z = (solid line) and z = 10£o (dashed line) at 
T = 0.2T c0 and 0.6T c0 where D = 20£ is fixed. 



the Fermi liquid corrections, 



Hf(k,r) 



= H 



h u + FS 



h u - 



1 M o 

3(1 + F a )Ff 
"2(3 + ffKt>pJVp 



M M (r) 
M N 

J* (r)L- (34) 



For realistic situation with Fq < 0, the enhancement of 
the surface magnetization M z (0) > Mn increases the ef- 
fective magnetic field, Hf ! (k,r)>H, while the suppres- 
sion of the magnetization in the middle region leads to 
Hf(k,r)<H. 

The field-dependence of the local spin susceptibil- 
ity M z (z)/Mn at finite temperatures is summarized in 
Fig. EJb). This ratio corresponds to the local spin sus- 
ceptibilities Xzz{z) for an applied field H \\f„, 



Xn M n 



(35) 



The spin susceptibility in the normal 3 He is obtained 
from xn = M N /H = 2^iV F /(l + F a ). In the low tem- 
perature regime, such as T = 0.2T c0 , the spin susceptibil- 
ities at the surface, Xzz(0)/xn, is enhanced in the linear 



regime of the magnetic field. As H increases, however, 
it reduces to xn at the surface, while Xzz{D/2) is insen- 
sitive to H as a result of the first-order transition from 
the B- to A- (or planar) phase. Although the nonlinear 
effect of the magnetic field suppresses the enhancement 
of the spin susceptibility due to the SABS at low tem- 
peratures, as seen in Fig. E^b), the spin susceptibility at 
the surface is still distinct from that in the central region 
of the system. In the higher temperature region where 
the B-phase undergoes the second-order transition to the 
planar or A-phase, Xzz(0) decreases to xn and Xzz(D/2) 
gradually increases as H increases. 



B. Relation between surface bound states and 
enhancement of magnetization 

As shown in Eq. ([TBI , a Zeeman magnetic field per- 
pendicular to the surface opens a finite energy gap. As 
shown in Fig[5j it simultaneously induces a large amount 
of the magnetization at the surface. Here, we clarify the 
relation between gapped surface bound states and the 
enhancement of the magnetization density. 

First, in Figs. Hta)-(c), we display the fc-resolved sur- 
face density of states, 



2 71 



N{k h r,E) = -hm I ^gK(k,r;E), 



(36) 



when the magnetic field is absent. In Eq. (|36|) . ku de- 
notes the momentum parallel to the surface, fc|| = sin0fc, 
as shown in Fig. [T] The retarded Green's function 
<7oHfe, r ; E) is obtained from Eq. (f2"2"j) with ioj n — > E+i0 + . 
Throughout this paper, we fix 0+ = 0.0057rT c o. Since 
the squashed B-phase in a slab geometry is SO(2)l + s 
symmetric around the i-axis, Af(k\\, r, E) describes the 
dispersion relation of the surface bound state. 

In the absence of a Zeeman magnetic field, the time- 
reversal symmetry as well as the particle-hole is pre- 
served. Hence, the BdG Hamiltonian H(k) is anti- 
commutablc with the chiral operator T combined with 
the time-reversal operator T = ia y T K and particle- 
hole operations C = cr x T y K, which is called the chi- 
ral symmetry, r} = 0. Here, K is the complex 
conjugate operator. The chiral symmetry allows one 
to introduce a three-dimensional winding number, w = 
I 2i^r e^Tr [L(1T 1 d^H) (H~ 1 d„U) (H~ 1 d v ffl\ , which is 
evaluated as w — 2 for the B-phase . 60 ' 62 Hence, the B- 
phase in the absence of a magnetic field is a topologi- 
cal phase and the bulk-edge correspondence implies the 
SABS satisfies E(k^) = at k\\ = 0, which is consistent 
with the analytic solution of the BdG equation within 
the Andreev approximation. 

The fc-resolved surface density of states for D = 20£o 
without a magnetic field, which is displayed in Fig.[6la), 
is consistent with the topological consideration, where 
the gapless point exists at feii = 0. However, since the 
SABS is localized at the surface within the coherence 
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length scale £o, the wavefunctions at both two surfaces 
are overlapped with each other as the thickness D ap- 
proaches £o- As discussed in Ref. I32l |79l - l82ll88l . the hy- 
bridization of wavefunctions localized at z = and D 
split the gapless cone as e~ D ^. Indeed, as seen in 
Figs-Ejb) andEJc), the spectral weight at fc|| =0 weak- 
ens as the thickness D approaches D C ri(0) = 9.6£o- In 
addition, it has the double peak in the low energy re- 
gion, where the upper branch has a distinct energy gap 
at fc|| =0 and another one remains almost linear at fi- 
nite fc||. For D = 10£o, the upper branch which has a 
energy gap E = 0.2ttT c q originates from the hybridization 
of Majorana cones bound at two surfaces, while the lower 
branch reflects the fact that the pair potential Aj_ which 
is perpendicular to the surface is squashed by two spec- 
ular surfaces as displayed in Fig. U) At D = D cr j(0), the 
squashed B-phase order parameter continuously turns to 
the planar phase with Aj_ = where fc|| =0 corresponds 
to the location of the point nodes in the bulk. The planar 
phase, the point node of which is normal to the surface, 
is not accompanied by the surface bound state and the 
low energy spectrum is linear on k\\ in the whole system. 

As seen in Fig. |Hl[d), the perpendicular field opens a 
finite energy gap in the surface cone, min \E\ ^0.157rT c o. 
For fi n H = 0.0488ttT c0 and T = 0.2T c0 , it is seen 
in Fig. [5] that the value of M z {z)/M^ at the surface 
z = is about 1.4. Then, the effective Zeeman en- 
ergy at the surface z = is estimated from Eq. (|3"4")l 
as Li n Hf(k,r)/(l + F a ) 0.2nT c0 . At z = 10£ , how- 
ever, it decreases to 0.l7rT c o, because of the suppression 
of the spin susceptibility M z (z = \Q^q)/M^k,QA. Hence, 
the energy gap min \E\ ^0.157rT c o in Fig.[5{d) is approxi- 
mately consistent with the analytic dispersion in Eq. (|16[) 
with the spatially averaged effective Zeeman energy. In 
the high magnetic field (Fig. EJe)), however, the nonlin- 
ear effect of the Zeeman magnetic field causes the pair 
breaking effect as displayed in Fig. [Jib). Therefore, as H 
increases, the bulk excitation gap becomes lower in ad- 
dition to the increase of the energy gap of surface bound 
state. This behavior is confirmed in Fig. [BJc) where 
the continuous excitation band lowers and merges to the 
gapped SABS branch. Since the situations of D = 12£o 
and 10£o in Figs. [6jf ) and (6jg) are close to the second- 
order phase transition field, the surface cone ceases to 
exist and the quasiparticle excitations in the entire sys- 
tem become gapless. 

Now, let us clarify how the change of the spectrum of 
the surface bound states affects the local magnetization 
density at the surface. In Figs. |7(a), [7£b), |7(d), we plot 
the surface density of states J\f(z = 0, E) in the absence 
and presence of a perpendicular magnetic field, respec- 
tively, which correspond to Fig. (Hla), HJd), and HJe). 
The local density of states Af(r, E) is defined as 

Af(r,E) = ~hn.($(k,r;E)). 

7T \ Ik 

= \\ d6 k sm8 k N(k\\,r,E), (37) 



(a) (b) (c) 
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FIG. 6: (Color online) fc-resolved surface density of states 
Af(k {[ ,z = 0,E) for £> = 20£o (a), 12£ (b), and 10£o (c), where 
H = 0. Af(k\\,z = 0,E) for D = 20Co at /i n # = 0.04887rT cO (d) 
and 0.08 547rT c0 (e). (f) and (g) are for D=12£ and 10£o at 
fi^H = 0.04887rT cO . In all the data, the temperature is set to 
be T = 0.2T c0 . 



where (•••)& denotes the average on the three- 
dimensional Fermi surface. The surface density of states 
in the absence of a magnetic field displayed in Fig. [TJa) 
is linear on E in the low energy region. The linear de- 
pendence reflects the gapless cone E(k)ot ^Jk^. + and 

is distinguishable from the full gap behavior at z = 10£o- 
In the case of H \\ z, Af(z = 0, E) is accompanied by the 
finite energy gap within \E\ <0.157rT c o. In the presence 
of a perpendicular Zeeman field, as seen in Figs.[7Jb) and 
[TJd), the finite energy gap appears in the low E region. 
In the high field regime, the surface density of states loses 
the linearity in the low E region and is indistinguishable 
from Af(z = D/2, E), due to the distortion of the order 
parameter induced by the nonlinear Zeeman effect. 

For comparison, we present in Fig. [SJa) the surface 
density of states with the gapless dispersion of the SABS 
under a magnetic field. This is realized when the mag- 
netic field is parallel to the surface {H \\ x) and n is fixed 
to be normal to the surface, n \ \ z, which corresponds to 
l z = R xz (n = z, ip) = in Eq. (fTo) . Hence, the surface 
density of states for H \\ x remains linear on E. Note 
that the configuration of ft \\ z becomes energetically un- 
stable in the strong magnetic field regime^ because the 
ground state has the n-vector texture which satisfies the 
condition £ z (n,(p) = R xz (n 7 tp) = l. The ground state un- 
der a strong parallel field is necessarily accompanied by 
the Zeeman energy gap of the surface bound state. 

Then, we introduce lm(g^(k, z; E))^, which is associ- 
ated with the contribution of quasiparticles in the super- 
fluid state to the local magnetization density M^z), that 



10 
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FIG. 7: Local density of states N(z, E) at z = (the solid line) 
and z = 10£ (the dashed line) for^ n # = (a), 0.0488ttT c o (b), 
and 0.08 547tT c o (d), where H is perpendicular to the surface. 
(c,e) — lm(g^(k, z; E))^ in the same condition as (b,d). In all 
the data, the temperature and thickness are set to be T = 
0.2T c0 and D = 20£ - 
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0.2 0.4 0.6 0.8 1.0 

(b) E/wTco 

FIG. 8: Local density of states N(z, E) at z — (the solid 
line) and 2 = 10{o (the dashed line) for ^ n # = 0.04887rT c o (b), 
where H || x and h \\ z. (b) — lm(g^ (fc, z; E))^ for the same 
situation as (a). The other parameters are same as those in 
Fig-El 

the bulk. 

This is in contrast to the case of the parallel magnetic 
field where the surface Majorana cone is assumed to re- 
main gapless. It is demonstrated in Fig. |SJb) that the 
low energy quasiparticles in the gapless Majorana cone, 
l-El < OAttTco, do not contribute to the magnetization. 
The contribution of — Im(g^-(fc, z; E)}^ in the higher en- 
ergy region E >0.6ttT c q comes up to the negative value, 
which suppresses the magnetization density relative to 
A/n- Hence, as long as the fi- vector is polarized to the 
surface normal, the magnetization density at the surface 
becomes highly anisotropic, which is associated with the 
dispersion of the surface bound state. This consequence 
is consistent with the interpretation of the surface bound 
states as the Majorana Ising spin. 



is, 

M„(z) - Mn^ =-^% i dElm (g*{k, z; -E)V . (38) 

fJL n H J \ Ik 

As seen from Fig. [TJc) with the solid line, the quantity 
— Im(<7^ (fc, z; E))^. at the surface z ~ becomes posi- 
tive in the energy region lower than the bulk excitation 
gap, E <0.67rT c o. This implies that the gapped surface 
bound state considerably enhances the local magnetiza- 
tion density, compared with Mn. As the Zeeman mag- 
netic field increases, however, the gapped surface bound 
state merges to the bulk excitations as seen in Fig. [SJe) 
and the positive contribution to M z decreases. This is 
seen in Fig. Etc) with the solid line. In the nonlinear 
regime of H, the resulting magnetization density at the 
surface becomes comparable with that in the normal 3 He. 
Note that since even within this regime the magnetiza- 
tion in the central region of the system (z = D/2) stays 
around M z ~ 0.4A'%, the enhancement of the magneti- 
zation at the surface is still distinguishable from that in 



V. PHASE DIAGRAM AND SPATIALLY 
AVERAGED SPIN SUSCEPTIBILITY 

In this section, we present the superfluid phase dia- 
gram of 3 He in a restricted geometry and the H- and 
T-dependencies of spin susceptibility averaged over the 
slab, where the latter is associated with the NMR fre- 
quency shift and absorption. In order to discuss the ther- 
modynamic stability and the phase diagram of superfluid 
3 He in a slab geometry, we estimate the thermodynamic 
functional within the quasiclassical approximation, 

^[5] = ^ 1 ^S P '^5 A -^)}, (39) 

where we set 

Sp'{---}=iV F J dr(Tr 4 {---}) W (40) 

The quasiclassical auxiliary function g is obtained from 
the quasiclassical Eilenberger equation (|22[) with replac- 



11 



ing £— >• \v_ (As [0, 1]), where the equation is solved once 
under a given self-energy but not self-consistcntly. The 
functional in Eq. (|39p is obtained from the Luttingcr- 
Ward thermodynamic functional associated with the 
Nambu-Gor'kov Green's function G, the detailed deriva- 
tion of which is followed by the work in Ref . l30l Equation 
(|39[) includes the influence of the condensation energy 
and quasiparticlc excitations as well as the Fermi liquid 
corrections. 



A. Effect of Fermi liquid corrections 

First, we emphasize that the thermodynamics is sen- 
sitive to the Fermi liquid corrections. Among their cor- 
rections, the Fq term associated with the local magneti- 
zation Mfj_(r) plays a crucial role. Figure [Hta) shows the 
field dependence of A^(z = D/2) of the squashed B-phase 
at Z? = 40£o~3.2/^m. It is seen from Fig. [9ta) that when 
the Fermi liquid corrections are absent, A±(z = D/2) con- 
tinuously vanishes at the critical field fi n H = 0.0827rT c o, 
where the second-order phase transition from B- to pla- 
nar (or A-) phase occurs. Figure [H£b) with the dashed 
line depicts the field-dcpcndcncc of the thermodynamic 
potential introduced in Eq. ([3"9"|) . where 5Hab(H,T) de- 
notes the thermodynamic potential of the B-phase rela- 
tive to the A-phase. 

As seen in Fig. IHIa) with solid lines, the Fermi liquid 
correction makes the field-dependence of A^(r) insensi- 
tive. This is because for a large D, the Fermi liquid cor- 
rection associated with M z (r) through Fft suppresses the 
effective magnetic Zeeman energy Hf(r)<H, except for 
the vicinity of the surface, as described in Eq. (fMf . As 
a result of the suppression of the effective field, the B- 
phase survives even in the higher magnetic field so that 
the phase transition from the B- to A-phase turns to the 
first-order transition, as displayed in Fig. Elb) with the 
solid line. 

As D decreases, however, the pair breaking effect at 
the specular surface gives rise to the distortion of the 
isotropic B-phase order parameter A± < A|| even in low 
fields and the spin susceptibilities become comparable to 
the value in the normal 3 He. In this case, the effective 
Zeeman energy /j, n H is unchanged from that of the bare 
Zeeman field fJ, n H and the Fermi liquid correction does 
not alter the qualitative feature of the phase transition. 



B. Phase diagram 

Figure QJJ summarizes the field and thickness depen- 
dences of Aj_(z = D/2) (solid lines) and the phase dia- 
gram (the bottom) at T = 0.2T c0 . In the region of the 
large thickness D > ll£o, the phase boundary is the first- 
order phase transition Hab- As D/£o increases, the first- 
order transition field Hab slightly increases and reaches 
saturation (J- n H AB /iTT c o = 0.095 in the thermodynamic 
limit fl>(o. Using the parameters T c q = lmK and the 
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FIG. 9: Field-dependence of A M (z)/7rT c0 at z = 20£ (a) and 
of thermodynamic potential of the B-phase relative to the A- 
phase, 8Q.AB (b), where D = 40£o and T = 0.2T c0 . The solid 
(dashed) line denotes SQ.\b with (without) the Fermi liquid 
corrections and the arrows in (b) point the first- and second- 
order transition fields. 
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FIG. 10: (Color online) Field-dependence of A±(z = D/2) for 
various values of the thickness D at T = 0.2ivT c o (solid lines). 
The open (filled) circles denote the first-order (second-order) 
phase transition points. The bottom describes the phase 
boundary between the distorted B-phase (the shaded area) 
and the planar (or A-) phase, where the thick (thin) line cor- 
responds to the first-order (second-order) line. 



gyromagnetic ratio of 3 He nuclei 7 = 2// n , the critical 
field is estimated as iJa R sa0.35T, which is consistent to 
Ref. HH and experiments in Rcfs. andl9ll. 

As the thickness D decreases the first-order transition 
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FIG. 11: (Color online) Superfluid phase diagram in the space 
spanned by temperature T, perpendicular magnetic field H, 
and thickness D. The definition of the open (filled) circles and 
thick (thin) lines are same as those in Fig. 1101 The shaded 
area is occupied by the distorted B-phase and the other is 
covered by the planar (or A-) phase. 



turns to the second-order, where it is seen from Fig. [10] 
that A±(z = D/2) continuously touches zero at the criti- 
cal field Hab- However, it is noted that the second-order 
phase transition may be proper to the weak coupling the- 
ory where the planar phase and A-phase are energetically 
degenerate. The finite contribution of an anisotropic in- 
teraction which makes A-phase more stable than the pla- 
nar phase, such as the spin- fluctuation feedback effect, 
may change the phase boundary to the first-order tran- 
sition even for small ZTs. 

In Fig. II H we summarize the phase diagram in a three- 
dimensional space spanned by the temperature T, per- 
pendicular magnetic field H, and thickness D. The first- 
order transition appears in the low temperature and large 
thickness region, while the high temperature and small 
thickness region involves the second-order phase transi- 
tion. Note again that the thermodynamic limit of this 
phase diagram, corresponding to D/^q^-oo, reproduces 
the well-known phase diagram in the bulk 3 Hc which 
is composed of the first-order transition in low T's and 
second-order line in high T's . 89 ' 90 

The bottom line in Fig. [11] describes the A-B phase 
transition in the absence of a magnetic field. The whole 
line is found to be the second-order transition, which re- 
produces the earlier works done by Hara and Nagai in 
Ref. and Vorontsov and Sauls in Ref. [33- However, 
Vorontsov and Sauls stated in Ref. |H that the second- 
order phase boundary around D ~ 10£o is covered by 
the new quantum crystalline phase, the so-called stripe 
phase, in which the translational symmetry in the plane 
of the film is spontaneously broken. In this paper, for 
simplicity, we eliminate the possibility of the stripe phase 
from the phase diagram. Since the stability against a 
magnetic field is not trivial, the complete phase diagram 
which takes account of the stripe phase remains as a fu- 
ture problem. 
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FIG. 12: (Color online) (a) T-dependence of the spatially 
averaged spin susceptibility (xzz) at fi n H/irT c o = 0.003 (solid 
lines), 0.03 (dashed lines), and 0.06 (dotted-dashed lines) for 
D = 12fo) 20£o, and 40£o- The arrows denote T m i n in the 
lowest field (see the text), (b), (c) T-dependence of local 
spin susceptibilities Xzz(z) at z = and D/2; (b) D = 20£o 
and (c) 12£ at fi n H = 0.0037rT c0 . The thin solid line in (a) 
and dashed line in (b) and (c) depicts the spin susceptibility 
X»* u m the bulk B-phase given in Eq. (|42l) . 



C. Spatially averaged spin susceptibilities 

Figure [Ufa) shows the T-dependence of spin suscepti- 
bility (xzz) averaged over the slab for D/£q = 12, 20, and 
40, where (xzz) is defined as 



1 

D 



D 



Xzz(z)dz. 



(41) 



For comparison, we plot the spin susceptibility in the 
bulk B-phase given with the Fermi liquid parameter Fq 
by 



v (buik) = (1 + F$)[2 + Y(T)] 
Xzz ' 3 + F a [2 + y(T)] AN ' 



(42) 



where Y(T) is the Yosida function^ The nonlinear effect 

of the Zecman magnetic field on y£ u was investigated 
by Fishman and Sauls in Ref. l93l 

According to the sum rule, the static spin susceptibil- 
ity (Xzz) is obtained by integrating the absorptive part of 
the dynamical spin susceptibility over all the frequency^ 
Hence, the temperature- and field-dependences are de- 
tectable through NMR experiments^ 

It is seen from Fig. IT2T a) that in weak magnetic fields, 
e.g., fi n H/irT C Q = 0.003, the spatially averaged spin sus- 



13 



ccptibility (xzz) has a minimum value at a certain tem- 
perature, e.g., T = T m i n «0.4T c0 for £> = 20£o, and raises 
up with further decreasing the temperature. This is con- 
trast to the T-dependence of the bulk B-phase associ- 
ated with the Yosida function. The increase of (Xzz) in 
the low-T regime is found to reflect the considerable con- 
tribution of the gapped surface bound state. To under- 
stand the non-monotonic behavior, we plot in Figs. 112( b) 
and [TWc) the T-depcndcnce of local spin susceptibili- 
ties Xzz{ z ) at z = and D/2, where we set D = 20£o 
in Fig. n^b) and 12£ in Fig. rT^c). In the case of 
D = 20£o, the spin susceptibility at z = 10£o traces the 
T-dependence of xfi u which corresponds to the bulk B- 
phase, and stays almost constant around Xzz{z~ D /2)~ 
0.4xn in low temperatures within T <0.4T c o. As we have 
discussed in Sec. IV B, a Zeeman magnetic field perpen- 
dicular to the surface opens the finite energy gap in the 
SABS, which gives rise to the enhancement of the spin 
susceptibility. In actual, it is seen from Fig. [T2T b) that 
Xzz(z) at the surface exceeds the Pauli susceptibility of 
the normal 3 He in the low temperature regime, while it 
monotonically decreases as T increases. Hence, the spin 
susceptibility (xzz) averaged over the slab in the low T 
region of Fig. [Ufa) indicates the enhancement of local 
magnetization density at the surface, while the behavior 
in the high T regime is dominated by the magnetization 
density in the central region of the system. 

As shown in Fig. 112( a). the qualitative feature on the 
T-dcpcndence of (xzz) is insensitive to the thickness D, 
except for the vicinity of the A-B transition Z?ab~9.6£o 
in which the magnetic response becomes indistinguish- 
able from that in the normal 3 He. It is demonstrated in 
Fig- [Ufa) with the solid lines that the temperature T m j n , 
at which (xzz) becomes minimum, lowers as D increases, 
namely, the thermodynamic limit is approached. 

Then, let us look at the field-dependence of (xzz), 
which is plotted in Fig.[T27a) with the dashed and dotted- 
dashed lines. In the bulk B-phase, as discussed in Refs.l87l 
andl93l. the nonlinear effect of the Zeeman magnetic field 
enhances the spin susceptibility in the entire region of 
T<T c0 as Xzz(H) -x Z2 (0)«(/ini?/A o (T)) 2 ^. The in- 
dependence also appears in the high T regime T > T m ; n 
of Fig. [l2T a). In the regime of T < T m ; n , however, the 
nonlinear effect of H lowers the magnetization density 
at the surface as shown in Fig. [5fb) and the resulting 
(Xzz) /xn is rather suppressed by increasing H. This im- 
plies that as the Zeeman magnetic field is ramped up, 
the characteristic temperature T m j n gets lower and fades 
away at last. In summary, the field and temperature de- 
pendences of (Xzz)/xn niay unveil the dispersion of the 
surface bound state in supcrfluid 3 He-B. 



VI. CONCLUDING REMARKS 

Here, we have investigated the role of surface bound 
states on the thermodynamics and spin susceptibilities 
in 3 He-B under a perpendicular magnetic field. First, 



within the Andrcev approximation of the Bogoliubov- 
de Gennes equation, we have clarified the relation be- 
tween the SO (3) order parameter manifold and the sur- 
face bound state, where the direction of the Majorana 
Ising spin is clarified. We have also explicitly mentioned 
that the condition in which the surface bound state opens 
the maximum energy gap coincides with the condition 
that makes the magnetic field energy lower. 

Subsequently, we have revealed the thermodynamics 
and surface bound states in a restricted geometry. All 
the results arc obtained with the quasiclassical Eilen- 
berger theory which provides the closed set of selfcon- 
sistent equations reliable to the weak coupling regime of 
superfluid 3 He. It turns out that the pair breaking effect 
and surface bound states play a crucial role on determin- 
ing the phase diagram and spin susceptibilities. The Zee- 
man magnetic field perpendicular to the surface always 
opens a finite energy gap in the surface bound state. We 
have demonstrated that the gapped surface bound state 
gives rise to the positive contribution to the enhancement 
of the spin susceptibility at the surface, compared with 
that in the normal 3 He. We have also emphasized the 
role of the Fermi liquid corrections in the phase diagram 
of a restricted geometry, which plays a critical role on 
determining the phase boundaries. 

We have also discussed the temperature- and field- 
dependences of the spatially averaged spin susceptibil- 
ity. It is found that the local spin susceptibility in the 
central region of the sample obeys the ordinary Yosida 
function, while at the surface it considerably increases 
in the low temperature regime. Hence, the temperature- 
dependence of the spatially averaged spin susceptibility 
in the low temperature regime is dominated by the con- 
tribution of the surface bound state, leading to the non- 
monotonic behavior. The characteristic temperature at 
which the spin susceptibility becomes minimum is sen- 
sitive to the thickness of the sample and monotonically 
decreases as the thickness increases. Furthermore, we 
have demonstrated that the nonlinear effect of the Zee- 
man magnetic field reduces the the spin susceptibility at 
the surface, resulting in the monotonic behavior on the 
temperature in the high field regime comparable with the 
A-B transition field. The spatially averaged spin suscep- 
tibility is detectable through NMR experiments. 

Finally, we would like to mention the issues of which 
we do not take account here: The effect of the surface 
boundary conditional— and the possibility of the stripe 
phase with the spontaneous breaking of the translational 
symmetry— In the absence of a magnetic field, the sur- 
face density of states in the low energy region is consid- 
erably enhanced by the diffusive surface. 30 ! 97 The low- 
energy density of states filled in by the skew scattering 
of the quasiparticle at the rough surface might drasti- 
cally change the temperature- and field-dependences of 
the spin susceptibility. Note that the specularity of the 
surface of 3 Hc can be experimentally controlled by coat- 
ing it with 4 He layers. 73 ' 74 Furthermore, the vicinity of 
the A-B phase transition around D ~ 10£ is occupied 
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by the stripe phased when the magnetic field is absent. 
However, the robustness against a Zeeman field is not 
trivial, which remains as a future problem. 
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Appendix A: Derivation of the dispersion in Eq. (|16[) 

Here, we describe the details about how to solve the 
Andreev equation (JSJ, 



-iav F k z d z r z + V + A( k F,a) <t>a{z) = Eip a (z), (Al) 



where <p±{z) describes the slowly varying part 
of quasiparticlc wavcfunction ¥>(?"), that is. 
f(r) = J2c^C a <f a (z)e lk ^ r with fc F , a = 
k F (cos^fc sin#fc , sin^fc sin6*fc , a cos k 9 k ) . The normal- 
ization condition is imposed on ip a (z) as 

Y, I <pl(z)<p a (z)dz = 1. (A2) 

ORf' 

The rigid boundary condition at 2 = 0, tp(x, y, z = 0) = 0, 
leads to C+ = — C_ and the continuity condition (p+(z = 
0) = ^_(z = 0). 

Then, introducing L£(n,ip) = diag[J7(n, <p), U*(n, <p)] 
and using the relation in Eq. ([?} , the BdG equation (|A1[) 
within the Andreev approximation reduces to 

[-iav F , z d z T z +V'(n,(f) + A (fc FjQ )] 

xU\h, <p)<p a (z) = EU\h, <p)<p a (z), (A3) 

where v F , z = v F cos8 k . The SU(2) matrix U_(n, <p) in the 
Nambu representation rotates the Pauli matrices and 
the Zeeman term V£' =lfi (n,<p)Vli(h,<p) results in 



V'(n,ip) = -finHfiRpvinffp) ( °" 



-a* 



(A4) 



It is convenient to introduce the unitary matrix M = 
(a x + a z )e i *< T */y/2 with Then, the pair 

potential Ao(fe) = ia li <TyAn,k f j, in the B-phase rotates to 
A^(fc F , Q ) = AfA (fc Fi± )M T , 



A' (k F , a ) 



where o,{4>hi^k) 



b{4>k,0k) -a*(<f>h,0k) 



(A5) 



AySin 2 4> k )sm6 k and b(4> k ,6 k ) 



A z cos6*fc + i(A x cos 2 (f> k 



-(A, 



that A x = A y = A z = Aq £ R. This is valid for a weak 
field regime within [i n H <C A , because the distortion 
induced by the magnetic field can be estimate in the 
thermodynamic limit as A||/A_l = 1 - O (/j, 2 H 2 / A 2 ) ^ 
where A|| (Aj_) represents the pair potential parallel 
(perpendicular) to an applied field. Within the as- 
sumption, Eq. (|A5[) reduces to a((f) k , 6 k ) — aAge k and 
b((t>k,9k) = 0. 

First, we solve the Andreev equation (|A1[) in the ab- 
sence of a magnetic field, H = 0. Then, the Andreev 
equation (|A3[) can be separated to two independent spin 
sectors as 



<Pa(z) 



fi 2 \z) 



&\z) 



(A6) 



where a/ and a g are the normalization constants. The 



equation for the wavefunctions f a 
Eq. CM) as 



(1.2) 



is obtained from 



Ho{k Fta ,z) 



(i) 

Di 

(2) 



En 



(1) 

Di 

(2) 



where 

U Q {k F , a ,z) = -iav FtZ d z a z + a x aA n e~ laSk ^ . 
Using the particle-hole symmetry, 

a x Ho{k, z)a x = -Ho(-k,z), 



(A7) 

(A8) 
(A9) 



the positive energy states with the wavefunction 
[fa,fai an( l E > is associated with the negative 
branch with ga' 2 ^ = <J x [fa , fa* 1 *] 1 " and —E. 

The resulting equation (|A7[) is equivalent to the one- 
dimensional Dirac equations with the mass domain wall. 
The index theore m 77 i"i 100 ensures the existence of the 
zero energy states when the mass term changes its sign. 
The bound state solution with ^(fen)! < Ao has the en- 
ergy dispersion linear on the momentum kn = (k x ,k y ) 
as 



£ (fc||) = ±^|fc|| 



(A10) 



This expression is independent of the orientation of n 
and the angle ip. The corresponding wavefunctions for 
the quasiparticles bound at at z = are given by 



/(fcj.,z)«(n > v)*±(0k) ) ( An ) 



where N k is the normalization constant estimated from 
Eq. (|A2|) . In Eq. (|A11|) , we also set f(k±,z) = 

sin (k±z) e~ z /£ with k± = . k F — kf, and 



-,±i- 



A y ) sin0fc cos(f> k s'm9 k . For simplicity, let us assume 





f 














1 




e 2 








% 






-i) 







(A12) 
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In Eq. (jAllI) . V'ofe corresponds to the positive energy 

solution and Vo~fe| * s the negative branch. The gapless 
spectrum of the SABS is protected by the nontrivial topo- 
logical invariant defined in the bulk region of the B-phase 
in the absence of a magnetic field.— ^ 

Now let us consider the case of a finite magnetic field 



tp k (r) = a+ip^l (r) + a-tpU (r), 



(A13) 



where the normalization condition for (r) requires 
\a + \ 2 + |a_| 2 = 1. The coefficients a± and energy E(k\\) 
are determined by solving the eigenvalue equation 



e^ fc 7z -\Eq\ 



a + 



I a a +_ ) , (A14) 



where 7 Z = fi n H^R^ z (n, <p) denotes the gap of the sur- 
face cone. From Eq. (|A14|) . the dispersion of the SABS 
is given as 



E(ku) = ±\ LEo(feii) + HnHt g (n,tp) 



(A15) 



and the wave functions are obtained from Eq. (|A13[) with 
a±(fc||) = 4/1(1 ± Here, we introduce the £- 

vector in Eq. (|A15|) . the definitio n 70 ' 75 of which is 



H 



(A16) 



Appendix B: Boundary conditions and numerical 
procedures 

The quasiclassical Green's function g is parameterized 
with 2x2 matrices a= a(fc, r; iu> n ) and b=b(k,r;iuj n ) as 



g(k, r; iu) n ) = -inN ( ^jt,^ 



2a 
-Co + ^ a 



where 



(0-0 - a&r 1 

{-ao + ba)- 1 



(Bl) 



(B2) 



This Ricatti parametrization automatically satisfies the 
normalization condition of g and simplifies the Eilen- 
berger equation (|22[) where the equations which gov- 
erns a and b are separated to each other. The resulting 
equations, called the matrix Ricatti equationsjiSir-iS^ are 
given by 

iv-p(k) ■ Va + 2iiu n a + A — aA^a + av' — va = 0(B3a) 
iv F (k) ■ Vb - 2iu n b + A f - bAb + bD- vb = 0,(B3b) 



where we set At = At(— k, r) and v = i/(k, r) is composed 
of the Fermi liquid correction and the Zeeman energy 



1 



v = v a + v^a^ 



l + F, 



1/ = v' Q a + i/> p - 



1 + Ft 



(B4a) 



(B4b) 



It is worth mentioning that the Ricatti amplitudes 
a and b have the following symmetry a(fc,r;icj„) = 
&*(— k, r; iuj n ), which implies that the quasiclassical 
Green's functions obey 



g j (k,r;i(j n )= g](-k,r;iu n ) 



f li {k,r;iw n )= /t( 



fc. 



r:iuj. 



(B5a) 



(B5b) 



where j — 0, x, y, z. 

In this work, we consider supcrfluid 3 Hc sandwiched 
by two specular surfaces, as displayed in Fig. [TJ As- 
suming spatial uniformity in the plane parallel to the 
surfaces, the resulting Ricatti equations (|B3|) reduce to 
one-dimensional ordinary differential equations along the 
z axis, which are numerically stable and requires an ini- 
tial value of a and b. For a quasiparticle momentum 
k on three-dimensional Fermi sphere with Vp(k) = v F k, 
we solve the Ricatti equations by numerically integrating 
along the classical forward (backward) trajectories for a 
(b) with an arbitrary initial value. The numerical inte- 
gration of one-dimensional Ricatti equations (|B3|) is per- 
formed with the fourth-order Runge-Kutta method from 
an arbitrary point of z. An arbitrary initial value of a 
and b converges after multiple reflections on the specular 
surfaces situated at z = and D. The specular surface 
requires the matching of two propagators a(k, z; ioj n ) and 
a(k, z; iu) n ), that is, 

a{k, z; iui n ) = a(k, z; iu> n ), for z = and D, (B6) 

where k = (cos cj)k sin#fc,sin^>fc sin#fc, — cos 9k). The 
boundary condition on b is given in the same way as 
Eq. (IB61) . 



Appendix C: SO (2) rotational invariance of the 
quasiclassical Green's functions 

Here, we clarify the symmetric property of the quasi- 
classical Green's functions. First of all, in the absence of 
a Zeeman magnetic field v = 0, the order parameter in the 
B-phase is isotropic in the sense of A M = Ao. Here, we 
introduce the simultaneous rotation in the orbital space 
Rflv = Ofj, u and spin space Rffl = (ROR" 1 )^ where 
{R)nv = Rp,v^ Using the rotation matrices, the momen- 
tum k and the Pauli matrices (or equivalently the 
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d-vector) behave as the three-dimensional vectors, which 
arc transformed to k^t-tR^ k v and h- >■ R fa? a v . Now 

let Us be an SU(2) representation of the 5*0(3) rotation 
(s) 

matrix . Then, the isotropic B-phase order parame- 
ter with A M = Ao is invariant under the joint rotation of 
spin and orbital spaces, SO(3)l+Si 

A(fe,r) = ia^yAoR^ku = U s A(R {L) k, r)E7j. (CI) 

However, a magnetic field and surface boundary con- 
dition reduces the joint rotational symmetry S0(3)l+s- 
Since we consider the situation where a magnetic field 
is applied along the surface normal (H _L z) and the 
dipole interaction is absent, it is natural to suppose that 
the components of the B-phase pair amplitudes still re- 
mains isotropic about the z-axis, A x = A y = Au and 
A z = A±_ =^ Am. Let O 1 - 2 ' be a two-dimensional rotation 
matrix around the z-axis and Ug be an SU (2) repre- 
sentation of the 50(2) rotation matrix {RO^RT 1 )^. 
Then, it turns out that the squashed B-phase order pa- 
rameter is invariant under the two-dimensional rotation 
of spin and orbital spaces around the z-axis, SO(2)l z+ s, j 



A(fc,r) 



^ 2) A(0 (2) fe,r)[7, 



(2)T 

s 



(C2) 



Now let us apply the SO(2)l x +s, rota- 
tion to the quasiclassical Green's functions as 
Z4aff (0^fc) f] iu) n )U2* = g(0^k,r\i(j n ), where we 

introduce U2 = diag[t/g , U^*] in the Nambu represen- 
tation. Then, the Pauli matrices and the momentum 
kfj, arc transformed to = [RO^ 1 R~ 1 )^ y a ll and 

k^ = 0$k u . The quasiclassical self-energy 5 in Eq. (|22|) 
is transformed to 1^S(0 {2) k,r)U^ = S(0^k,z) in the 
same way. The magnetic Zeeman term in Eq. (1221) with 
H = Hz is invariant under the SO(2)l^ + s z rotation. 
To this end, the Eilenberger equation (|2"2"]l under the 
S0(2)l z +s z rotation reduces to 

j n r z -S{0 {2) k 1 z)-v,g{0 {2) k 1 z- ) iuj n ) 

+iv F k z d z g{0 {2) k, z; iu n ) = 0. (C3) 

Under the SO(2)l :c+ s z rotation, the quasiclassical self- 
energy matrix S(k, z) mapped to 



HO$k,z) A(k,z) 
At(-fc,z) tf{-0$k,z) 



, (C4) 



where we use Eq. (|C2|) and the diagonal part is 

£(0 (2) fc, z) = ^ (O (2) fe, z)a + ^(0 (2) fe, z)?„. (C5) 



As described in Eq. (|25[) . the terms vq and are ex- 
panded in terms of the Legendre polynomials Pg. Among 
the possible contributions, we suppose in this paper that 
only the £ = and 1 channels play a crucial role on ther- 
modynamics and surface bound states, and the contri- 
butions with the higher £'s are eliminated. For 3 He-B, 
one finds 1^0 = because of the absence of the mass flow. 
Hence, Eq. (|C5|) reduces to 

v(pV*k,z) = -fi^lH^-M^z)}^ 

Now, we suppose that the quasiclassical self-energy v is 
invariant under the SO(2)l z +s z rotation in 3 He-B with 
a perpendicular magnetic field, 



£>(0 (2) fe,z) = u(k,z). 



(C7) 



This requires that the magnetization density and the 
spin current J must satisfy the following conditions: 

M x (z) = M y (z) = 0, J xy {z) = -J yx (z). (C8) 

With the sclfconsistent calculation of the quasiclassical 
Eilenberger equations, we confirmed that the B-phase 
under a perpendicular magnetic field always satisfies the 
conditions in Eq. (|C8|) . 

To this end, the Eilenberger equation within the 
SO(2)l,+s z symmetry is written as 



iuj n r z - S(k, z) - w,g(0 (2) fc, z; iu n ) 

+iv F k z d z g(0 {2) k, z; iu n ) = 0, (C9) 



which gives the equation for the quasiclassical Green's 
function at a point 0^ 2 'k of the three-dimensional Fermi 
sphere. Equation (|C9[) is also equivalent to the equation 
for g(k, z; iu> n ). Namely, there is a one-to-one correspon- 
dence of the quasiclassical Green's functions between two 
points k and Ok and the quasiclassical Green's function 
at an arbitrary point of the Fermi sphere, O^k, is ob- 
tained from g(k, z; iu3 n ) as 



g(0^k, z; ioj n ) = W 2 g{k, z; iuj n )U_ 2 



(CIO) 



This relation through the 50(2) rotation is useful for 
shorting the computation time of the selfconsistent cal- 
culation. Once we calculate g(k, z;ioj n ) along the path 
(i) displayed in Fig. [3j the Green's function g for all 

k is given by the symmetric relation in (|C10[) with 
g(k,z;iuj n ). 
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